library(sf)
library(ggplot2)
library(dplyr)
library(tmap)

Mekansal Analiz nedir?

Mekansal analiz icin gerekli olan materyaller nelerdir?

Shape File ve Ozellikleri

-Mekansal analiz için en temel gereksinimlerden biri shape file denilen geometri dosyasıdır. - Basit anlamda shape file geometrik bilgiyi ve o geometriye ilişkin coğrafi özellikleri içerir. -Bir shape file’daki geometrik bilgiler nokta, çizgi ya da poligon ile temsil edilebilirler. -Uygulamada basitçe bir ülkedeki şehirlerin her birinin birer poligon olarak ifade edilebilir.

# turkey.shp ==> Shape File
# turkey.shx ==> Autodesk AutoCAD tarafindan derlenen, sekil ve font tanimlamalarini iceren dosya
# turkey.dbf ==> data base file

# prj, 
#.gda, 
#.sbx, vb.

Turkiye’nin Shape File’i

https://data.humdata.org/dataset/turkey-administrative-boundaries-levels-0-1-2

Shape File Yukleme

turkey <- st_read("C:/Users/firat/Desktop/turkey_administrativelevels0_1_2/adm1a/tur_polbnda_adm1a.shp")
## Reading layer `tur_polbnda_adm1a' from data source `C:\Users\firat\Desktop\turkey_administrativelevels0_1_2\adm1a\tur_polbnda_adm1a.shp' using driver `ESRI Shapefile'
## Simple feature collection with 81 features and 8 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: 25.66851 ymin: 35.80842 xmax: 44.81793 ymax: 42.10479
## Geodetic CRS:  WGS 84
str(turkey)
## Classes 'sf' and 'data.frame':   81 obs. of  9 variables:
##  $ adm1_tr   : chr  "ADANA" "ADIYAMAN" "AFYONKARAHISAR" "AGRI" ...
##  $ adm1_en   : chr  "ADANA" "ADIYAMAN" "AFYONKARAHISAR" "AGRI" ...
##  $ adm1      : chr  "TUR001" "TUR002" "TUR003" "TUR004" ...
##  $ Shape_Leng: num  8.27 5.49 7.15 6.59 4.43 ...
##  $ Shape_Area: num  1.407 0.751 1.45 1.163 0.79 ...
##  $ adm0_en   : chr  "TURKEY" "TURKEY" "TURKEY" "TURKEY" ...
##  $ adm0_tr   : chr  "TÜRKIYE" "TÜRKIYE" "TÜRKIYE" "TÜRKIYE" ...
##  $ adm0      : chr  "TUR" "TUR" "TUR" "TUR" ...
##  $ geometry  :sfc_MULTIPOLYGON of length 81; first list element: List of 1
##   ..$ :List of 1
##   .. ..$ : num [1:2863, 1:2] 36.3 36.3 36.3 36.3 36.3 ...
##   ..- attr(*, "class")= chr [1:3] "XY" "MULTIPOLYGON" "sfg"
##  - attr(*, "sf_column")= chr "geometry"
##  - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA NA NA NA NA NA NA
##   ..- attr(*, "names")= chr [1:8] "adm1_tr" "adm1_en" "adm1" "Shape_Leng" ...
plot(turkey)

plot(st_geometry(turkey))

Harita Projeksiyonu

st_crs(turkey)
## Coordinate Reference System:
##   User input: WGS 84 
##   wkt:
## GEOGCRS["WGS 84",
##     DATUM["World Geodetic System 1984",
##         ELLIPSOID["WGS 84",6378137,298.257223563,
##             LENGTHUNIT["metre",1]]],
##     PRIMEM["Greenwich",0,
##         ANGLEUNIT["degree",0.0174532925199433]],
##     CS[ellipsoidal,2],
##         AXIS["latitude",north,
##             ORDER[1],
##             ANGLEUNIT["degree",0.0174532925199433]],
##         AXIS["longitude",east,
##             ORDER[2],
##             ANGLEUNIT["degree",0.0174532925199433]],
##     ID["EPSG",4326]]

Burada son satırda gördüğümüz “ID["EPSG",4326]]” bilgisi haritamızın projekte edilmiş olduğunu gösteriyor. Peki “st_crs” kodunu çalıştırdıktan sonra şöyle bir ifade ile karşılaşmış olsaydık ne yapacaktık?

#proj4string: "+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs

Bu ifade bize haritamızın projekte edilmediğini gösterecekti çünkü “proj4string” ifadesini “+proj=longlat” takip ediyor. Eğer karşılaştığımız sonuç bu şekilde olsaydı, o zaman haritamızı projekte etmemiz gerekecekti. Peki hangi projeksiyon bilgisini kullanmamız lazım? Bu sorunun cevabi derin felsefi tartışmalar gerektirir ki bizim amacımız sadece haritamızı doğru bir şekilde görselleştirmek. Bunun için ülkelerin EPSG kodlarını bulmamız isimizi görecektir. Bunu yapabilmek için de ülkelerin UTM paftalarını bilmemiz gerekiyor.

Kısa bir Google araştırmasından sonra Türkiye’nin 4 farklı UTM paftasına sahip olduğunu görüyoruz.

Turkiye’nin UTM Paftasi